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PRODUCTION  AND  CONTROL  OF  ION  CYCLOTRON  INSTABILITIES 
IN  THE  HIGH  LATITUDE  IONOSPHERE  BY  HIGH  POWER 
RADIO  WAVES 

I.  INTRCDUCTICN 

It  is  now  recognized  that  many  of  the  phenomena  occurring  in  the 
auroral  ionosphere  may  oe  associated  with  the  presence  of  the  field-aligned 
current-systems  which  are  a  constant  feature  at  these  latitudes  [see  e.g., 
Burke,  1981].  In  addition  to  being  a  key  component  of  the  ionosphere- 
magnetosphere  coupling  circuit,  these  currents  may  also  be  responsible  in 
generating  micro-and/or  macro-structures  (turbulence)  in  the  ambient 
plasmas  at  low  altitudes.  Such  structures  have  been  observed  by  various 
measurement  techniques  (ground  based  ionosonde,  VHF  radar  backscatter, 
optical  measurements,  scintillations,  in-situ  rocket  and  satellite 

measurements ,  etc.)  A  catalogue  of  the  structure  scale  sizes  spans  the 
domain  from  a  few  meters  to  a  few  hundreds  of  kilometers  [see  e.g.,  Aarons, 
1973;  Clark  and  Raitt,  1976;  Dyson,  1969;  Dyson  and  Winningham,  197*1; 
Fremouw  et  al.,  1977;  Greenwald,  1974;  Hanuise  et  al.,  1981;  Bower  et  al., 
1966;  Kelley  et  al.,  1980;  Ogawa  et  al.,  1976;.  Olesen  et  al.,  1976;  Phelps 
and  Sagaiyr.,  1976;  Sagalyn  et  al.,  1974;  Vickrey  et  al.,  1980;  Weaver, 
’965].  Recent  theoretical  studies  have  attempted  to  interpret  these 

structures  in  terms  of  plasma  instabilities,  neutral  turbulence,  structured 
particle  precipitation,  convection  flows,  etc.  [Xeskinen  and  Ossakow, 
198;].  Due  to  the  ever-varying  conditions  prevalent  in  the  auroral 

ionosphere,  it  seems  likely  that  all  of  these  processes  may  be  occur i.ng  at 
one  time  or  another.  In  some  instances,  different  processes  can  result  in 
similar  observable  effects  which  makes  the  resolution  of  individual 

mechanisms  difficult  and  introduces  complications  in  the  modelling 

process.  The  use  of  active  experiments  in  the  r.ear-earth  space  offers  an 
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attractive  complement  to  the  studies  of  understanding  of  natural  plasma 

processes  occurring  in  the  ionosphere.  One  technique  that  has  been  widely 
used  is  the  modification  of  tne  ionosphere  using  powerful  radio  waves 

[Utlaut,  1970;  Fialer,  1974;  Thide  et  al . ,  1982;  Stubbe  et  al . ,  1981;  Wong 
et  al . ,  1931]  A  considerable  amount  of  theoretical  work  has  been  done  in 
this  field  [Perkins  and  Kaw,  1971;  Perkins  et  al . ,  1974;  Fejer,  1979; 

Vaskov  and  Gurevich,  1977;  Das  and  Fejer,  1979;  Stenflo,  1983].  One  of  the 
possible  consequences  of  the  interaction  of  a  large  amplitude  EM  wave  with 
the  ionosphere  is  the  generation  of  structures  in  the  ambient  plasma 
density.  This  may  come  about  several  ways;  one  of  them  being  the 
excitation  of  plasma  instabilities  by  an  electromagnetic  pump  wave  under 
sub-threshold  conditions.  By  the  same  token,  one  could  use  the  same 
principle  of  parametric  coupling  of  the  external  pump  and  the  natural 
plasma  modes  to  suppress  a  current-dr iven  or  gradient-dr iven  plasma 
instability  in  the  system.  Work  on  basic  parametric  instability  processes 
in  plasmas  is  exhaustive  [Silin,  1965;  DuBois  and  Goldman,  1965;  Nishikawa, 
1968].  The  application  of  these  ideas  to  the  possible  control  of 
ionospheric  irregularities  was  originally  made  by  Lee  et  al.  [1972]  for  the 
case  of  equatorial  electrojet;  and,  has  since  been  applied  for  the 
equatorial  spread  F  situation  [Bujarbarua  and  Sen,  1973]  and  the  auroral 
ionosphere  [Kesxinen  et  al.,  1983].  In  this  report,  we  examine  the 
possibility  of  external  control  (excitation  or  suppression)  of  the  current- 
driven  ion-cyclotron  instability  in  the  auroral  ionosphere.  In  addition  to 
being  excited  in  the  topside  ionosphere  [Kindel  and  Kennel,  1971],  the 
current-driven  ion  cyclotron  instability  may  also  be  triggered  in  the 
collisions!  bottomside  ionosphere  [Chaturvedi,  1976;  Satyanar ayana  et  al., 
1435;  Fejer  et  al.,  1984]  by  a  radio-wave  at  the  local  upper-hybrid 
frequency.  In  the  next  section,  we  outline  the  theory,  and,  estimate  the 


J 
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incident  power  density  required  to  stabilize  or  destabilize  the 


instability.  The  final  section  considers  application  of  these  results  to 
the  bottomside  high  latitude  F  region  ionosphere. 


II.  THEORY 

We  follow  the  theoretical  procedure  outlined  in  Lee  et  al.  [1972]  and 
Xeskinen  et  al.  [1983]  for  the  calculation  of  the  dynamic  plasma  response 
to  a  large  amplitude  electromagnetic  pump  wave.  The  equilibrium  consists 
of  a  field-aligned  current  Jqz(=  "  n0eVQ“)  and  an  oscillating  electric 
field  E0  =  E  cos  unt.  We  have  equated  the  zero-order  current  to  an 
equilibrium  electron  drift  V and  have  used  the  so-called  dipole 
approximation  for  the  pump  wave  which  is  valid  when  the  wavelength  of  the 
perturbations  ( —  k  ^  )  is  considered  much  smaller  than  the  scalesize  of  the 
variation  of  the  pumo-field  (k^1),  i.e.,  kQ  <<  k.  We  shall  use  the  two- 
fluid  equations  in  describing  the  electron  -  and  ion  -  fluid  dynamics, 
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velocity  of  lig.nt,  electric  field  and  magnetic  field,  respectively.  Since 

we  restrict  this  study  to  bottom-side  altitudes  we  take  v. .  <  v. 

i  i  -  i  n 

where  v.  .  is  the  ion-ion  collision  frequency  and  v.  the  ion-neutral 
collision  frequency.  In  equilibrium,  we  have 
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where  a  =  e,i,  and  _7q  and  refer  to  drifts  induced  by  the  ambient  and 

pump  fields,  respectively.  Note  that  eq.  (7)  describes  the  approximate 

equilibrium  oscillatory  motion  of  plasma  particles  under  the  action  of  the 

?.F  electromagnetic  field.  Linearizing  (1)  -  (4)  for  perturbations, 
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It  is  convenient  to  transform  to  an  oscillating  frame,  defined  by 


The  plasma  quantities  in  the  oscillating  frame  are 

'r'  ,  tl  =  of  (r  ,  t]  where  the  orime  denotes  the  oscillating 
a '  2  '  a  at  ' 

quantities.  Equations  (3)  -  (9)  are  similar  in  the  oscillating  frame, 

substitutions  fv  .  +  V  )  +  V  5n  on  ,  ov  -*•  ov  ,  5E  ♦  5E  and  r  -*■ 
aO  oi  p  au  a  oi  —a  —a  —  —  — 

Assuming  a  harmonic  perturbation,  exp  [i(k  •  r^  -  mt) ]  with  _k  =  (k^x  * 
k,  y  *  k^z)  and  »  =  u>r  *  iY,  we  may  combine  (3)  and  (9)  and  write 
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We  have  used  =  0,  7.  =  0  in  deriving  (13)  ~  ( 1 a .  One  now  transforms 

0  l 

back  to  the  laboratory  frame  and  writes  down  the  dispersion  relation  for 
the  modes  by  using  the  Poisson’s  equation.  This  transformation  has  been 
discussed  in  bee  et  al.  I ‘972],  and  the  dispersion  relation,  thus  obtained, 
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ana  ;  is  defined  by  k  •  3  =  ;  sin^t  +  3 )  ,  and,  Hg  =  -  -j  (l  +  4w 

hot  a  tnat  we  have  ignored  the  effect  of  the  pump  wave  on  ions  and 
have  use-i  trie  approximation  ;  <<  1  in  deriving  (15)  which  3tates  that  the 
electron  excursion  length  in  the  pump  field  is  smaller  than  the 
perturbation  wavelength.  In  absence  of  the  pump  field,  t.ne  left  hand  siae 
of  ,’5.  leads  ta  the  dispersion  relation  for  ion-cyclotron  waves 
CChaturvedi ,  * 97b I , 
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where 


=  T  /m;  .  An  application  of  this  instability  for  tne  auroral 


Ionosphere  has  recently  been  discussed  by  Satyanarayana  et  al.  [1935]. 


Physically,  the  instability  is  related  to  the  dissipative 
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electrons  in  their  parallel  motion  which  impedes  the  "instantaneous" 
redistribution  of  the  electron  fluid  to  a  Boltzmann-like  distribution  in 
the  wave  potential.  In  the  collisionless  limit,  the  dissipation  is  due  to 
t.ne  wave-particle  (landau)  resonance  [Drummond  and  Rosenbluth,  1962],  while 
in  the  collisional  case,  it  is  caused  by  the  electron  collisions  with  ions 
'  or  neutrals)  [Chaturvedi  ana  Kaw,  1975].  When  the  electron  drift  velocity 
exceeds  the  parallel  wave  phase  velocity,  the  electrostatic  ion  cyclotron 
(SIC)  wave  becomes  unstable. 

In  the  presence  of  the  pump-wave,  the  dispersion  relation  (17)  is 
modified  to 
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used  2  /v  >  k/k  ,to  >  Y,  v.  ,  and,  '  >  A/fiT,  in  deriving  '20''  -  .'21  ).  We 
can  rewrite  the  modified  growth  rate  of  the  mcdes  as 
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where  Y'  denotes  the  linear  growth  contribution  in  absence  of  the  pump- 
field  while  Y,  represents  the  modification  introduced  by  the  pumo.  We  have 

3 

again  used  ;  <<  1  in  going  from  (21)  to  (25).  We  see  from  inspection  of 
'25)  that  the  presence  of  the  pump  can  be  stabilizing  (destabilizing) 


5  >  0  ( <  C),  or,  when  <  w,,,,].  We  note  that  similar  criteria 
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were  also  obtained  by  earlier  workers  [Lee  et  al .  ,  1 972,  Sujarbarua  and 

Sen,  1978,  Xeskinen  et  al.,  1983]. 

We  can  treat  the  finite  ion  temperature  effects  in  a  similar 
fashion.  Since  the  ion-cyclotron  instability  has  maximum  growth  rate  in 


the  ionosphere  near  -  1,  [Satyanarayana  et  al . ,  1935;  Xindel  and 


Kennel,  1971]  one  needs  to  use  the  kinetic  (collisionless  Ylasov-Boltzmann) 
equation  for  ions  [Dun  and  Dupree,  1970].  In  this  case,  the  modified 
dispersion  relation  becomes. 
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:ne  icse.noe  of  tr.e  pump  wave,  this  expression  yields  the  ion-cyciotrcn 
wave  distension  relation. 
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yields  the  expressions  for  the  real  frequency  and  the  growth  rati 
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where  the  second  term  in  (29)  in  the  expression  for  Y  results  from  the  ion- 
lar.dau  damping  and  we  have  added  a  contribution  (last  term)  due  to  ion- 
viscous  damping  [Oum  and  Dupree,  1970]  for  completeness.  The  collisions! 
ion-cyclotron  instability  described  by  (29)  may  have  relevance  for  the 
short  wavelength  observations  (such  that  k ,  p ^  -  1  )  of  ion-cyclotron  waves 
at  low  ionospheric  altitudes  [Bering,  198d;  Satyanarayana  et  al.,  1935]. 

The  inclusion  of  pump  wave  effects  in  equation  (23)  leads  to 
modification  of  (29)  in  the  form  given  by 
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where  A  and  Y.  are  given  by  (22)  and  (29)  respectively.  We  see  that  the 
pomo  wave  is  stabilizing  (destabilizing)  for  5  >  0  (<  0),  i.e., 
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One  can 


when  a)  >  id. <  'd,.,.>;  as  in  case  of  fluid  ions  (eq.  (25)). 
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readily  extend  this  analysis  to  the  case  of  the  kinetic  EIC  instability 
[Drummond  and  Rosenbluth,  1962].  However,  in  the  ionosphere,  the  kinetic 
instability  is  likely  to  occur  at  higher  altitudes  [Kindel  ana  Kennel , 
1971]  such  as  above  the  F  peak;  whereas,  the  modification  experiments  have 
access  only  to  the  lower  collisional  alti.udes. 

A  rough  estimate  for  the  threshold  amplitude  of  the  pump-field 
required  to  stabilize  or  destabilize  the  collisional  EIC  instability, 
reveals  that  (using  (30)), 
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where  t  is  given  by  (27)  and  we  have  considered  a  circularly-polarized  0- 
ode,  so  ip  =  -  m/2  and  j  |  =  |E„„  j.  For  typical  high-latitude  F-region 


py  1 


parameters  [Satyanarayana  et  al.,  1985],  we  have,  3^  -  0.5  G, 


-9 


va/P.a  ~  10  >  iTJ  ~  10  S  '  1  >5 ,  pL  -  1 0^cm ,  id0  -  wpo  -  3  x 

T  =7.  =0.1  eV,  p  =  3  cm,  V  -  2  x  10J  cm  s  1 ,  v.  =0.1  sec  1 , 

3  1  Q  J  1 


x  10 


7  -1 


v.  .  =  0.31  sec  ,  giving  t.nat 


-J 


'  „  _  -4  2 

/Sir  =o  x  10  W/m  . 


III.  DISCUSSION 

We  have  studied  the  effects  of  a  high  frequency  electromagnetic  pump 
wave  on  the  evolution  of  ion  cyclotron  instabilities  in  a  weakly  ionized 
collisions!  plasma.  We  find  that  an  electromagnetic  pump  wave  oscillating 


with  the  upper  hybrid  frequency  can  lead  to  the  artiricial  production  or 

control  of  ion  cyclotron  instabilities  in  a  ccllisional  plasma.  The  basic 

physical  mechanism  of  this  process  can  be  expressed  as  follows.  The 

electromagnetic  pump  wave  with  frequency  interacts  with  low  frequency 

ion  cvclotron  modes  co  and  drives  sideoands  at  a  +  u>MU.  These  sidebands 

—  Jn 

can,  in  turn ,  couple  to  the  pump  fields  and  nonlinearly  affect  the  ion 

cyclotron  modes.  This  interaction  is  ponderomotive  in  nature.  However, 

the  partial  pressure  force  (Fejer,  1979)  can  also  play  a  role  depending 

upon  the  k  of  the  sidebands.  Fejer  [1979]  has  shown  that  the  partial 

pressure  force  exceeds  the  ponderomotive  force  roughly 

when  a  ">  -U  where  X  is  the  electron  mean  free  oath.  For  the 
\  e*  e 

collisions!  SIC,  (k_/k  )  -  0.06,  ( k^p . )  -  1.5  and  -  1  km.  [Satyanarayana 

et  al.,  1985],  which  implies,  X  <  X a .  Thus,  for  the  ion-cyclotron  modes, 

the  oredominant  pump-wave-induced  parametric  effects  are  likely  to  be  of 

ponderomotive  force-like  in  nature.  For  parameters  typical  of  the  high 

latitude  ionosphere,  we  find  that  ion  cyclotron  modes  can  be  stabilized  or 

destabilized  with  an  incident  pump  electric  field  on  the  order  of  0.1 

7/m.  The  polarization  of  the  pump  wave  should  be  0  mode.  These  power 

density  levels  could  be  accessible  using  the  ionospheric  heaters  in  Alaska 

(Wong  et  al.,  1981)  and  Norway  (Stubbe  et  al.,  1981). 

We  have  made  several  simplifying  assumptions  in  the  course  of  this 

paper.  We  have  assumed  the  dipole  approximation,  i.e.,  taken  the 

wavelength  of  the  oumD  wave  -  k  1  to  be  much  larger  than  the  wavelength  of 

P 

the  ion  cyclotron  modes  k  1 .  Since  k  1  -  100's  m  and  k  -  p.  -10m,  this 

P  i 

assumption  is  valid.  In  addition,  we  have  neglected  pump  electric  field 
swelling  near  the  reflection  point  (Ginzburg,  196-0  and,  as  a  result,  our 


threshold  oumo  electric  fields  should  be  considered  as  lower  bounds. 
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